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Using the RAVE survey, we recently brought to light a gradient in the mean galacto- 
centric radial velocity of stars in the extended solar neighbourhood. This gradient likely 
originates from non-axisymmetric perturbations of the potential, among which a per- 
turbation by spiral arms is a possible explanation. Here, we apply the traditional dcn- 
Vh \ sity wave theory and analytically model the radial component of the two-dimensional 

velocity field. Provided that the radial velocity gradient is caused by relatively long- 
lived spiral arms that can affect stars substantially above the plane, this analytic 
model provides new independent estimates for the parameters of the Milky Way spi- 
ral structure. Our analysis favours a two-armed perturbation with the Sun close to the 
inner ultra-harmonic 4:1 resonance, with a pattern speed f2 p = 18.6lg'2 kms" 1 kpc^ 1 

and a small amplitude A = 0.55^o 02% °f ^ ne background potential (14% of the back- 
ground density). This model can serve as a basis for numerical simulations in three 
dimensions, additionally including a possible influence of the galactic bar and/or other 
non-axisymmetric modes. 

Key words: Stars: kinematics - Galaxy: fundamental parameters - Galaxy: kine- 
matics and dynamics. 



1 INTRODUCTION 

It has long been recognized that internal secular evolution 
processes should play a major role in shaping galaxy disks. 
Among the main drivers of this secular evolution are the 
* E-mail: arnaud.siebert@astro.unistra.fr disc instabilities and associated non-axisymmetric pertur- 
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bations, including the bar and spiral arms. Quest ions about 
their nature - transient or quasi-stationary (e.g. , ISellwoodl 
l2010al : iQuillen et al. IlioTll : brand et al. ll2012ah - , about 
their detailed structure and dynamics such as their am- 
plitude, pattern speed, pitch angle or number of arms, as 
well as questions about their de tailed influence on secular 
processes like stellar migr ation (jSellwood fc Binnev 1 I20021 ; 
iMinchev fc Famaev I I2OI0I) . are all essential elements for a 
better understanding of galactic evolution. The Milky Way 
provides a unique laboratory in which a snapshot of the dy- 
namical effect of present-day disc non-axisymmetries can be 
studied in great detail, and help answering the above ques- 
tions. 

Current knowledge of the structure and dynamics of 
the bar and the spiral arms of the Milky Way relies 
both on the gas, a nd notably on i ts ob se rved longitude- 
veloc it y diagram l|Birmev et al~l Il99ll ; iBissantz et al. I 
120031 ; Englmaier et al. l201ll) or maser s in high mass 
star forming regions dReid et al. I 120091. and references 
therein), and on the stars (e.g.. iGeorgelin fc Georgelin I 



19761: iBinnev. Gerhard fc Spergel II 19971; IStanek et al. Ill997l ; 
Benjamin et al. 20051 ; iLepine et al. Il2011bl ). For the spiral 
arms, bo th type s of co nstraints are combined in a recent 
study by IValleel (|2008l ). whose model predicts the location 
in space and velocity for the spiral arms. 

With the advent of new spectroscopic and astrometric 
surveys, six-dimensional phase-space information for stars 
in an increasingly large volume around the Sun allow us 
to set new dynamical constraints on the non-axisymmetric 
perturbations of the Galactic potential. An example of such 
new detailed kinematical information on stellar motions in 
the extended solar neighbourhood is the recently detected 
(ga lactocentric) radial vel ocity gradient of ~ 4kms~ 1 kpc -1 
by ISiebert et al~l <|2011al ), making use of more than 200 
thousand stars from the RAVE survey. If this result is not 
owing to systematic distance errors (which the geometry of 
the radial velocity flow seems to exclude by not depend- 
ing on distance and longitude in any systematically biased 
way), and more importantly, if one assumes that, at first 
order, what is seen above the plane is a reflection of what 
would happen in a razor-thin disc, and that the spiral arms 
are long-lived, one can apply the ana lytic density wave de- 
scription of spiral arms proposed by iLin fc Shu I (|1964| ) to 
constrain the shape, amplitude and dynamics of spiral arms. 

Whether long-lived density waves are the correct de- 
scription of spiral patterns in galaxies remains heavily de- 
bated. From a theoretical point of view, while it seems 
that the radial velocity dispersion profile nee ded to support 
long- lived spiral waves in bar less discs (e.g . . iBertin fc Linl 
1 19961 ) would be heavily unstable (|Sellwood Il2010al ). the sit- 
uation is much less clear in the presence of a central bar, 
where nonlinear mode coupling between the bar and spi- 



ral could sustain a long-lived spiral pattern JVoglis et al 
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20121 ). while iGrand et al. I (|2012bl ) however find that spiral 
arms are transie nt even in the pre s ence o f a central bar. On 
the other hand, iD'Onghia et al. 1 (|2012l ) find locally long- 
lived self-perpetuating spiral arms which could be locally 
consistent with density waves, but fluctuating in amplitude 
with time. Furthermore, long-lived spirals can also develop 
as being sustained by coherent oscilla tions due to a flyby 
galaxy encounter l|Struck et al. Il201ll ). a process we know 



to be ongoing for the Milky Way, and cosmologically simu- 
lated disk galaxies exhibit a distribution of young stars con- 
sistent with the pred ictions of classical dens ity wave theory 
for long-lived spirals (|Pilkington et al~ll2012l ). Finally, let us 
note that both long-lived and transient spirals can coexist 
in a galaxy, which adds complexity to the picture. 

On the observational side the situation is also un- 
clear. Evidence seems to exist for both transient and 
long-lived spiral arms: e.g., M81 apparently contains long- 
lived spiral arms consistent with the classical density 
wave theory jLowe et al. Ill994l ; lAdler fc Wefstpfahl I Il996l ; 
iKendall et al. 1120081 ) while in M51, even if its disc streaming 
motion appears consistent with the density wave des cription, 
the m ass fluxes are inconsistent with a steady flow (|Shettv I 
120071 ). Studying observational tracers for different stages of 
the star-formation sequence in 12 nearby spiral galaxies, 
iFovle et al. I |201ll ) also found that they do not show the 
expected spatial ordering for long-lived spiral arms, from up- 
stream to downstream in the corotating frame. In the Milky 
Way, many of the dynamical constraints on spiral arms cur- 
rently come from local constraints provided b y velocity space 
substructures also k nown as moving groups l|Dehnen Ill998l ; 
iFamaev et al~ll2005l ). For instance, examining the local stel- 
lar distribution in action space, ISellwood I l|2010bl ) found 
that stars from the Hyades moving group were concentrated 
along a resonance line in action space, which was interpreted 
as a signature of scattering at the inn er Lindblad reson ance 
of a transient spiral pattern (see also [McMillan l201ll . who 
found that this feature could also be associated with an 
outer Lindblad resonance). However, in this picture, only 
the Hyades moving group is accounted for, and the remain- 
ing substructures observed in the local phase space distri- 
bution must be explained by invoking oth er origins. Models 
based only on transient spiral arms (e.g.. |Pe Simone et al. I 
12004 ) were actually unable to reproduce the precise loca- 
tion of the various other prominent moving groups, such as 
Sirius. On the other hand, models based on long-lived spi- 
ral arms, locating the 4:1 inner resonance close to the Sun, 
were able to reproduce both the positio n of the Hyades and 
Siriu s moving groups at the same time (|Quillen fc Minchev I 
l2005l; iPompeia et al. II20TTI ) as well as other moving groups 



I Antoia et al. 2011). Other observational arguments based 



on the step-like metallicity gradient in the Galactic disc also 
argue in favor of long-lived spirals (|Lepine et al. Il201ll ). 

Given this theoretical and observational situation, we 
here make the conservative assumption that interesting in- 
formation can be retrieved from the classical analytic treat- 
ment of spiral arms as long-lived density waves. This ana- 
lytic model could then serve as a basis for numerical simu- 
lations in a three-dimensional disc. The paper is structured 
as follows. In Section 2, we review the data, as well as the 
analytic density wave model we are using. We present and 
discuss our results in Section 3, and conclude in Section 4. 



2 DATA & METHOD 

2.1 Two-dimensional velocity field 

Our analysis is based on data from the RAVE survey 
JSteinmetz et al. Il2004 IZwitter et al~ll200Sl ; ISiebert et atl 
2011b) which provides line-of-sight velocities with a preci- 
sion of 2 kms _1 for a large number of bright stars in the 
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southern hemisphere with 9 < I < 12. RAVE selects its tar- 
gets randomly in the I-band interval, and so its properties 
are similar to a magnitude limited survey. The RAVE cat- 
alogue is cross-matched with astrometric (PPMX, UCAC2, 
Tycho-2) and photometric catalogues (2MASS, DENIS) to 
provide additional proper motions and magnitudes. In this 
study, we use the internal version of the catalogue which 
contains data for 434,807 spectra (393,903 stars). To com- 
pute the galactocentric velocities, knowledge of the dis- 
tance to the star is required. For R AVE stars, distances t o 
30% are available in thre e studies: Breddels et al. I (2010), 
IZwitter et al~l (|2010h and iBurnett et al. I (|201lh ~All cata- 
logues provide compatible distance estimators and the ve- 
locity maps generated using the different catalogues are sim- 
ilar. 

Our final sample consists of 213,713 stars from this sur- 
vey limited to a distance of 2 kpc from the Sun and to a 
height of 1 kpc above and below the plane. We demon- 
strated the existence of a velocity gradient of disc stars 
in the fou rth quadrant, direct ed radially from the Galac- 
tic centre l|Siebert et al. Il2011al ). The two-dimensional mean 
galactocentric radial velocity field in the Galactic plane is 
presented in Fig. [T] where we use a box 4x4 kpc in size, 
centered on the Sun, sampled using 60 bins in each direc- 
tions. For this analysis, we restrict ourselves to bins con- 
taining more than 5 stars and the mean velocity is com- 
puted using a median function. Note that converting ve- 
locities in the heliocentric reference frame into the (Vr, Ve) 
galactocentric coordinates requires the galactocentric radius 
of the Sun _Ro, the Sun's peculiar velocity vq with respect 
to the Local Standard of Rest (LSR), and the motion of the 
LSR with respect to the Galactic center vlsr- We assume 
Rq = 8 kpc for the distance of the Sun to the Galactic cen- 
tre and Vlsr = 220 km s -1 to match the values used for the 
mass model (see Section 12. 2p . We use th e latest determi- 
nation of the value of the solar motion by ISchonrich et al. I 
(|201GT ): (Uq, Vq) = (11.1, 12.24) kms" 1 . 

The gradient affects a sample dominated at large dis- 
tances by red giants, with a typical velocity dispersion 
<jr ~ SO^Okms" 1 , and affects stars substantially above the 
plane, keeping in mind that RAVE lines of sight are typically 
at |b| >~ 20°. The zone where the gradient is the steepest 
is populated with stars with typically \z\ ~ 500 pc. How- 
ever, if ?7lsr is positive (a local mean motion towards the 
inner Galaxy), then (Vr) would be by construction negative 
in the Sun's neighbourhood and reach at larger distances 
and larger heights. In the modelling procedure hereafter, we 
let {Vr)r = — UhSR be a parameter of the model to allow 
us to account for uncertainties on this quantity. 

Ideally one would also use the tangential velocity field 
(Ve) in combination to the (Vr) field. However, as stated 
above, our sample reaches distances to the Galactic plane of 
1 kpc, avoiding the regions close to the plane. The RAVE 
survey mimicking a magnitude limited survey in fields 6 deg 
in diameter on the sky, each field containing a different num- 
ber of stars, the stellar population mixture varies from point 
to point on the maps. Therefore, the contribution of the 
asymmetric drift is difficult to estimate while it enters the 
calculation of the Ve component. Hence, we chose to restrict 
our analysis to the (Vr) field, although we give the full set 
of equations, including this component, in the next section. 



2.2 Density Wave model 

To model the velocity field, we use th e density w a ve de - 
scription of spiral arms p roposed by ILIn fc Shu | (Il964 



and further developed in Lin. Yuan fc Shu 1 19691 ) and 



IShu. Stachnik fc Yostl (]l97ll ). This model is based on 
an asymptotic analysis of the WKBJ type of the Eu- 
ler/Boltzmann equations, valid only in the regime of weak, 
long-lived and tightly wound spirals (small pitch angle). This 
model being well known a nd documented (see for example 
iBinnev fc Tremaine Il200"8h . we restrict its description to the 
main results used in this study. 

The perturbation to the potential considered is of the 
form 



$i = A(R) exp [i(ut - m0 + #(#))] 



(1) 



where A(R) is the amplitude of the perturbation, m is the 
number of arms and <&(R) is a monotonic function. The per- 
turbation rotates at an angular frequency given by 



■ u/m . 



(2) 



The perturbations in the components of the mean ve- 
locities (Vr) and (Ve), where (R,0) are the coordinates in 
the cylindrical coordinate system centered on the Galaxy, 
are given by 



(Vr) = 
(Ve) = 



kA v 

K 1 — v' 2 

1 kA 



2 n l-v 2 



(x) cos(x) 
(x) sin(x) 



where 



(3) 



(4) 



k being the radial wave number, or the velocity dispersion, k 
the epicyclic frequency and v is defined by v = m(Q p — Q) /k. 

The functions J-l and J-l are the "reduction factors" 
that correct the mean velocities for the effect of velocity 
dispersion, lowering the effect of the spiral perturbation on 
the velocity field as the velocity dispersion increases. In the 



limit of a zero velocity dispersion (J-J = T v — 1) we 
recover the velocity field of the gas while if or becomes large, 
Ju = — ¥ and the velocity field becomes unaffected 
by the spiral p erturbation. The two fun ctions are given in 
Appendix B of iLin. Yuan fc Shu I (| 19691 ) ■ 
The phase of the spiral pattern \ is defined by 



X = ut - m9 + $(R) 



(5) 



which, in the case of logarithmic spirals with &(R) = 
mcotgi lni? in Eq.Q](i being the pitch angle), can be writ- 
ten in a more convenient form 



X = Xo + m(cotgi ln(R/R ) -(6- 8 )) . 



(6) 



The subscripts in the previous equation refer to the value 
at the Sun's location. The radial wave number k is then 
given by 



k(R) = $'(#) = mcotgi/R, 



(7) 



with k(R) < for trailing waves and k(R) > for leading 
waves. 

The mean velocities of Eq.[3]vary in the Galactic plane 
as a function of R and 8, the modulation depending on the 
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Figure 1. Left: velocity field for the radial component of the velocity vector (Vr) as a function of location in the Galactic plane. Right: 
associated random error on the mean velocity. In both panels, the location of the Sun is marked by a star, the Galactic centre is towards 
positive x and the y-axis is oriented towards the Galactic rotation. The colour coding follows the mean velocity and mean velocity error 
in kms -1 . 



mass model via Q and k, the velocity dispersion in the radial 
direction or via the reduction factors and on the parameters 
of the spiral perturbation. These parameters are the number 
of arms m, the amplitude of the perturbation A, its pattern 
speed Q p , the pitch angle i and the phase Xo- In addition we 
chose to include (Vr)r as a free parameter while comput- 
ing the model to account for possible uncertainties on Ulsr- 
While this parameter is usually taken into account while 
computing the velocities (Vr, V$) of the observations, here 
we include it as a correction to the predicted (Vr) and (Vg) 
in the model to avoid the computation of the velocity field 
at each step which is time consumming. In the remainder of 
the paper, we will denote P the vector of model parameters 
P = (m,A,n p ,i,Xo,(TR, (Vr)r ). 

For the rotatio n curve of the Milky Way, we use the 
models I and II of iBinnev fc Tremainel teOO&t) Table 2.3 , 
based on the mass models of iDehnen fc Binnev I (|l998h . 
These models reproduce equally well the circular-speed 
curve and other observationally constrained quantities such 
as the Oort constants, the surface mass density within 1.1 
kpc or the total mass within 100 kpc of the Milky Way. The 
two models correspond to two limiting cases where either the 
disc or the halo dominates the rotation curve (model I and 
II respectively). The models being computed for Ro — 8 kpc 
and V co ~ 220 kms" 1 , we use the same values for these two 
parameters when computing the galactocentric velocities in 
Fig. ID 

In practice, our tests showed that all our solutions con- 
verged on approximately the same value for or. This is 
due to our sample being dominated by the old thin disc 
population and we chose to fix or to the best fit value, 
or = 34.2 km s -1 , to reduce the dimension of our parameter 
space. This value of the dispersion fits very well the observed 
dispersion from the RAVE sample, excluding the tails rep- 
resentative of the thick disc population and of large proper 



motion errors. Hence, the model parameters we consider in 
the remainder of the paper is P = (m, A, fi p , i, xo, (Vr)r ). 

The mean velocities of Eq. [3] are compared to the two 
dimensional velocity field of Section \2. II The comparison is 
done using a chi-square estimator 



2 _ ( (VR)i,obs - (Vr) i, model ) 



where the sum is on all bins containing at least 5 stars 
and <7j l0 f, s is the error on the mean velocity shown in Fig. [T] 
right panel. We restrict the analysis to the mean velocity in 
the radial direction, the tangential velocities being affected 
by the asymmetric drift which can not be properly taken into 
account in the model, our sample being a mixture of stellar 
populations of different ages, even though it is dominated 
by the old thin disc (see Section [271} . 

We note also that systematic distance errors would af- 
fect the results presen ted below. As shown in the first paper 
ijSiebert et al. ll2011al ). a systematic error in the distances 
affects the measured gradient in (Vr) by approximately the 
same factor: a 20% overestimate/underestimate of the dis- 
tances induces a ~ 20% overestimate/understimate of the 
amplitude of the velocity gradient in (Vr) , which to first 
order, results in a higher/lower amplitude A of the spiral 
perturbation by the same amount. However, as shown in 
the same paper, an independent estimate of the velocity 
field using red clump stars, for which an unbiased distance 
estimate can be obtained, shows a good agreement of the 
velocity fields, giving us confidence that our distances can 
not be strongly affected by an unknown bias and we will not 
consider this possibility in the remainder in this paper. 
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3 RESULTS AND DISCUSSION 
3.1 Parameter space sampling 

The number of arms in the Milky Way is not known with 
certainty. Both 2-armed and 4-armed spiral pattern are con- 
sidered in the literature although an m = 2 mode in the 
stars seems to be favoured. In the analysis we will consider 
both cases and look for the best matching solution for either 
number of spiral arms. 

Also, some recent works have suggested that the local 
standard of rest is not on a circular orbit, eg. Ulsb. might not 



be kms~ 1 (ISmith et al. 1120091: ISmith Whiteoak fe Evans I 
l2012l ; lMoni-Bidm. Carraro fe Mendez \\20lt l. Therefore, we 
include the possibility for a non-null U component of the 
LSR in the model. 

For the minimisation, we considered the standard 
value of [/lsr = Okms" 1 ^] as well as values of 
±5 kms~ 1 whose amplitud e corr espond to the finding of 
ISmith. Whiteoak fe Evansl (|2012r i. Finally we left {Vr)r 
as a free parameter in the fit. However, we note that our sam- 
ple reaches only 2 kpc away from the Sun, not deep enough 
in the plane to disentangle the effect of a radial motion of 
the LSR from uncertainties in the U component of the so- 
lar motion with respect to the LSR. Therefore we should 
keep in mind that a non-null best fit value of the {Vr)r 
parameter does not necessarily imply a radial motion of the 
LSR. 

The summary of the chi-square analysis is presented in 
Table Q] and the chi-square contours for the best models are 
presented in Fig. [2] In this figure, the two panels show the 
1, 2 and 3-a contours, fixing all the other parameters to 
the best fit solution, in the Q p versus amplitude plane (left 
panel) and pitch angle i versus phase xo (right panel). The 
plain lines are for the mass model I, the red dashed lines for 
model II. 

The best fit model is obtained for a two armed spi- 
ral mode with the mass model II. The best fit solution for 
model I is equally good with a chi-square difference of 0.5. 
Generally, the difference between the two mass models is 
low. This is expected as within the region sampled by the 
RAVE data, the models are comparable with fi m i — f2 m ii « 
0.35 kms" 1 kpc -1 . The density wave model being not sen- 
sitive to the details of the mass model - the latter entering 
the equations only through Q(R) and k— the two models can 
only be distinguished in regions where they are significantly 
different (for example closer to the Galactic centre). 

The chi-square value of the best m = 4 solution is also 
close to the best m = 2 solution. However, for m = 4 
the pitch angle i is found to be « —26 degrees, out of 
the bounds of the tight-winding approximation upon which 
the density wave model relies (max ~ 15 — 20 degrees, 
iLin. Yuan fc Shu I (| 19691 ) ). Hence, for a four-armed pattern, 
we conclude that no satisfactory solution is found and we 
will discard four armed patterns in the following discussion. 

Focusing on the m = 2 mode, a strong correlation 
is observed between the amplitude and the pattern speed 
while the correlation is weaker between the pitch angle and 
the phase (Fig. If the phase and pitch angle are well 
determined, the shape of the contours in amplitude ver- 



Recall that C/lsr : 



■ (Vr)r ■ 



sus Q p indicate a large range in possible solutions at the 
3-<t level: Q p varies from 12 to 22 kms" 1 kpc _1 and the 
amplitude from 0.1 to 0.9% of the background potential. 
The amplitude of the best fit model is A — 0.55% of the 
background potential. This value translates to 14% of the 
ba ckground density wh i ch is close to the value proposed 
by iMinchev fc Famaevl (|2010h and consistent with earlier 
determinations summarized in lAntoia et al. I (|201 lh for the 
local spiral amplitude. We also note that corotating waves 
({l p ~ 27.5 kms" 1 kpc" 1 ) seem to be excluded. 

Finally, among the m = 2 solutions, a low radial compo- 
nent of the LSR velocity is prefered. The best fit model con- 
verges to (Vr)r = 0.9kms _1 while a zero radial compo- 
nent can not be ruled out when comparing the chi-quare val- 
ues. On the other hand, a more pronounced outwards motion 
of the LSR ( (Vr)r = 5 kms" 1 ) shows significantly larger 
chi-square value s while an inwards motion of — 5 kms" 1 
as suggested by ISmith. Whiteoak fc Evans I (|2012l ) is even 
less consistent. However in the latter case, the model con- 
verges outside of the range of allowed values for the pitch 
angle which limits the conclusions one can draw from this 
resu lt. We note that our b est fit value is consistent with 
the ISchonrich et al. I (|2010l ) errors at the 2-cr level when 
considering the error bars on the determination of U (e.g. 
U — 11.1°'q 9 75 kms -1 ). In the next section, we will concen- 
trate on the best m = 2 models and their implications for 
the structure of the Galactic disc. 



3.2 Resonances and spiral structure 

The best fit (Vh) velocity fields for the m = 2 solutions 
using the mass models II is presented in Fig. [3] Although 
both reproduce equally well the structure of the velocity 
field between y = ±1 kpc, a small difference between the 
two solutions exists, model I predicting larger velocities in 
the top right and lower left corners of our sample (bottom 
right panel). However the velocity difference reaches only 
0.2 kms -1 , a level that lies below the capabilities of our 
data. The region below y = — 1 kpc and at < x < 1 kpc is 
apprently poorly reproduced, however the lower left panel 
of Fig. [3] indicates that our solution stays well within the 
observational errors. The right panel of Fig. [T] indicates that 
this region suffers from large velocity errors and has there- 
fore a lower weight in the solution. In the regions where our 
data are of the best quality (mostly \y\ < 1 kpc), our models 
reproduce adequately the observed velocity field. 

Focusing on the pattern speed, our best models sug- 
gest that the Sun is located ~ 200 pc inside the in- 
ner 4:1 resonance (ultra-harmonic renonance or UHR) of 
the spiral pattern (Fig. [4]). Our finding for the pattern 
speed Q p =18- 19km s" 1 kpc" 1 is in agreement with re- 
cent studies that also place the Sun clos e to th e UHR: 
Q p = 17 kms" 1 k pc" 1 by lAntoia et al. I dgOllT) . Q p = 
18 km s " 1 kpc" 1 bvlOuillen fc Minchev I l|2005h or Q p /Q = 
0.65 by iPompeia et al7 ( 201lh to be com pared to Q v /Qq ~ 
0.68 in our study. However, as shown by [Gerhard] l|201ll) , 
determinations of the spiral arms' pattern speed range from 
17 to 28 kms" 1 kpc" 1 , the higher values being prefered by 
open cluster birthplaces while hydrodynamical simulations 
and phase space substructures favour slower pattern speeds. 
It is interesting to note that the pattern speeds found from 
velocity space substructures (jQuillen fc Minchev I 120051 : 
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Table 1. Chi-square results. Parameters with error bars were kept free in the minimisation. The error bars correspond to the 1-cr 
internal errors obtained from the chi-square contour. Models marked with (*) have a large pitch angle (open arms) and do not satisfy 
the tight-winding approximation. The number of pixels used in the minimisation procedure is 1595. 
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Figure 2. Cuts through the chi-square space around the best fit models in the amplitude vs pattern speed plane (left panel) and phase 
vs pitch angle plane (right panel). The crosses mark the location of the best fit while the contours are the 1, 2 and 3 cr limits. The 
amplitude in the left panel is given as a percentage of the background potential at the Sun's location. The plain lines are for the mass 
model I and the dashed lines for the mass model II. In both panels the top contours are for m = 4, the bottom contours for m = 2. 



lAntoia et al. I l201ll ; iPompeia et ai~ll201ll ) are close to our 
value. This would indicate a similar origin for the velocity 
gradient and the velocity substructures, reinforcing our as- 
sumption that the velocity gradient we observed is due to 
spiral arms. However this statement must be put in perspec- 
tive as both types of study rely on the same assumptions 
that the spiral pattern is long-lived and tightly wound. 

Comparing the predicted den sity pattern to the loca - 
tion of the spiral arms obtained bv lEnglmaier et al. I (|201lh 
in the gas, we find a good agreement (Fig. U right panel). 
Both the Perseus arm and the Centaurus arm are recov- 
ered at the proper location. The Sagittarius-Carina arm is 
not recovered in our models. This indicates that this fea- 
ture is not a dominant feature in the Solar neighborhood, 
reinforcing the view that the Milky Way spiral arm pat- 



tern is dominated by two main arms, Perseus and Centaurus 
(iDrimmel ll2000l:lDrimmel fc SpergeT1l200ll ; lBeniamin et al. I 
l2005l ; IChurchwell et al. 1120091 ). Contrary to the m = 2 mod- 
els, our best fit m = 4 solution does not reproduce any of 
the known spiral arms, in addition to being invalidated by 
its large pitch angle, hence an m = 2 mode for the spiral 
pattern in the Milky Way is preferred within the Lin-Shu 
regime. 



4 CONCLUSION 

We have an a lysed the velocity gradient detected by 
ISiebert et al~l l|2011al ) using the R AVE data in the fr ame- 
work of the density wave model of iLin fc Shu I (| 19641 ). as- 
suming that the velocity gradient we detected is due only to 
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Model II 




x (kpc) x (kpc) 



(model-obs)/a Model I - Model II 




x (kpc) x (kpc) 

Figure 3. Top panels: observed velocity field (left) and model velocity field for the best fit solution using the mass model II (right). 
Bottom right: velocity field difference between the best fit solutions using the mass model I and II. The colour coding follows the 
median galactocentric radial velocity in kms" 1 in the three panels. Bottom left : difference {Vn m ll) ~ (Vr observed) normalised by the 
observational errors showing that all velocities on the 2 dimentional map are well recovered within the observational uncertainties. On 
all panels, the Sun's location is at (0,0) and is marked by the black or white star. 



spiral arms and that the spiral arms in the Milky Way are 
long-lived. 

Our model converges properly for an m = 2 pattern, 
while if the chi-square of the m = 4 solutions are compara- 
ble, the predicted pitch angle is too large, invalidating the 
solution. 

The best fit solutions for m — 2 reproduce adequately 
the observed velocity field for (Vr) in the region |y| < 1 kpc 
where our data are the most reliable. Outside of this region, 
the difference between the model and the observations is 
still within the observational errors although the agreement 
is less clear. 

The predicted pattern speed places the Sun about 
200 pc outside the inner UHR of the spiral arms. Such a lo- 
cation of the Sun is consistent with previous works based on 
velocity space substructures, suggesting a similar origin for 
the velocity space substructures and the (Vr) gradient. Our 
best fit value for the amplitude of the spiral perturbation, 
A = 0.55+q;o2% of the background potential or 14% of the 



backgro und density, is consi s tent w ith the value proposed 
by, e.g.. iMinchev fc Famaevl (|2010l ) and is al so in the range 
of ear lier measurements as summarized in lAntoia et al. I 
l|201ll ). 

Comparing our model to the location of spiral arms in 
the gas, we find a good agre ement with the lo c ation of the 
major spiral arms given by lEnglmaier et al.l l |201lf ). The 
density enhancement predicted by our best model matches 
the location of the Perseus and Centaurus arms. The Sagit- 
tarius arm is not reproduced by our solution which tends to 
reinforce previous studies concluding that the Milky Way 
spiral potential is dominated by a two-armed mode, the 
Sagittarius-Carina arm being a minor feature for the dy- 
namics of the disc. 

Our study r elies on the density wave model of 
iLin fc Shu I (ll964T ) and we assumed no vertical variation 
of the (Vr) field within the limit of our data. RAVE data 
do contain the 3-dimensional spatial information which 
we will use in further studies. However, going from 2D 
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Model I 



Model I 




Figure 4. Left: circular frequency as a function of galactocentric distance for the mass model I of Binney & Tremaine (thick lines). 
The dotted and dashed lines are respectively the relations Q ± re/4 and ± re/2 versus R. The location of the Sun is marked by the 
vertical dotted line at Rq = 8 kpc . The horizontal lines are the pattern speed corresponding to the best fit models for m = 2 (plain lines) 
and m = 4 (dash-dotted lines). Right: Density associated to the best model. The grey shading and contours represent the overdensity 
associated to the spiral perturbation. The contours are evenly spaced by O.lSo, the background column density, from 0.1 to 0,5£q. The 
blue c ontour depicts the footprint of the RAVE data. The red lines mark the location of the spiral arms in the gas from lEnglmaier et al. I 
(2011). From left to right we have Cygnus arm (top left corner), the Perseus arm, the Sagittarius-Carina arm and the Scutum-Centaurus 
arm. The results for the mass model II are almost identical and are therefore not presented. 



to 3D requires an upgrade of our modeling technique 
taking properly into account the asymmetric drift and 
the vertical variation of the spiral potential. Moreover, 
vertical p erturbations leading to possible variations of 
(V z )(R.z) dSrnith. Whiteoak fe Evansll2012l ; IWidrow et al. 1 
|2012| ; IWilliams et al. I lin prep]) are intrinsically not taken 
into account in our analysis. Future 3D simulations of such 
perturbations and their possible influence on the (Vr) field 
will be necessary to disentangle their possible effects from 
the velocity gradient modeled here. Finally we note that our 
model is local as a result of th e tight- winding approximation 
(see for example discussion in lBinne"v"ll2012l Section 1.4.2). 
Ongoing surveys like Gaia-ESO (GES) or SDSS/SEGUE 
will provide data in the Galactic plane that can be used 
to test our models further in towards the Galactic centre 
(GES) or further out (SDSS/SEGUE). It will be interesting 
to test whether these two surveys predict the same pattern 
speed for the spiral arms. 
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